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Abstract 

Let M be a smooth compact oriented Riemannian manifold of dimension n without boundary, 
and let A be the Laplace-Beltrami operator on M. Say / / 6 5(K + ), and that /(0) = 0. For 
t > 0, let K t (x,y) denote the kernel of f(t 2 A). Suppose / satisfies Daubechies' criterion, and b > 0. 
For each j, write M as a disjoint union of measurable sets Ej z k with diameter at most ba J , and 
measure comparable to (bet 3 )™ if ba? is sufficiently small. Take Xj t k £ Ej,k- We then show that the 
functions cj>j t k(x) = ^(Ej^k) 1 ^ 2 K aJ {xj t k, x) form a frame for (I — P)L 2 (M), for b sufficiently small 
(here P is the projection onto the constant functions). Moreover, we show that the ratio of the frame 
bounds approaches 1 nearly quadratically as the dilation parameter approaches 1, so that the frame 
quickly becomes nearly tight (for b sufficiently small). Moreover, based upon how well-localized a 
function F G (I — P)L 2 is in space and in frequency, we can describe which terms in the summation 
F ~ SF = ~^2 k (F, 4>j,k)4>j,k are so small that they can be neglected. If n = 2 and M is the torus 
or the sphere, and f(s) = se~ s (the "Mexican hat" situation), we obtain two explicit approximate 
formulas for the 4>j,k, one to be used when t is large, and one to be used when t is small. 
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1 Introduction 

Say /o € (the space of restrictions to E+ of functions in 5(M)). Say f ^ 0, and let 

/(*) = sfo(s). 

One then has the Calderon formula: if c e (0, oo) is defined by 

r°° Ht r 00 
c = J o \f(t)\ 2 j = J o t\f Q (t)\ 2 dt, 

then for all s > 0, 

r°° dt 

I \f(t S )\ 2 j = c<^. (1) 
The even function h(£) = /(£ 2 ) is then in <S(R) and satisfies 

(In fact, all even functions in S(M) satisfying (2) arise in this manner). 

Its inverse Fourier transform h is admissible (i.e. is a continuous wavelet). We prefer to write, formally, 

h = f(-d 2 /dx 2 )S; 

the formal justfication being that 

(f(-d 2 /dx 2 )sy = f(e) = h(o. 

Thus f(—d 2 /dx 2 )S is a continuous wavelet on R. 

Discretizing (1), if a ^ 1 is sufficiently close to 1, one obtains a special form of Daubechies' condition: for 
all s > 0, 

OO 

0<A a < \f(a 2j s)\ 2 <B a <cx., (3) 

j=-oo 

where 

Aa ^2]ToH (1 " ° (l(a " 1)2(1 ° S 10 " ' W 

Ba = ^oga\ ^ + ° il{a ~ 1)2(1 ° g |a _ • (5) 

((4) and (5) were proved in [8], Lemma 7.6, if a > 1. If a < 1, replace a by l/a in that lemma.) In 
particular, B a /A a converges nearly quadratically to 1 as a — > 1. For example, Daubechies calculated that 
if f(s) = se~ s and a = 2 1 / 3 , then B a /A a = 1.0000 to four significant digits. 

Calderon's formula and Daubechies' condition are very important in the construction of continuous 
wavelets and frames on BL 

Our program is to construct (nearly tight) frames, and analogues of continuous wavelets, on much more 
general spaces, by replacing the positive number s in (1) and (3) by a positive self-adjoint operator T on 
a Hilbert space H. If P is the projection onto the null space of T, by the spectral theorem we obtain the 
relations 
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r°° Ht 

j o \f\ 2 {tT)- = c{I-P) (6) 

and 

CO 

A a (I-P)< \f\ 2 {a^T)<B a {I-P). (7) 

j=-oo 

(The integral in (6) and the sum in (7) converge strongly. In (7), X^l-oo := ^ m M,N^oo SjL-M> taken 
in the strong operator topology.) (6) and (7) were justified in Section 2 of our earlier article [9]. (This 
article should be regarded as a sequel of [9], although, if one is willing to take some facts (including (6) 
and (7)) for granted, it can be read independently of [9].) 

Taking T to be —d 2 /dx 2 on R leads to the continuous wavelet f(—d 2 /dx 2 )S on R. (Of course, on R, 
P = 0.) 

We began our program of looking at more general positive self-adjoint operators T, in our article [8]. 
There we took T to be the sublaplacian L on L 2 (G), where G is a stratifed Lie group, and thereby 
obtained continuous wavelets and frames on such G. 

In this article we will look at the (much more practical!) situation in which T is the Laplace-Beltrami 
operator on L 2 (M), where (M, g) is a smooth compact oriented Riemannian manifold without boundary, 
of dimension n. We will construct nearly tight frames in this context. Now P will be the projection onto 
the one-dimensional space of constant functions. We constructed continuous wavelets on M in [9]. 
We discussed the history of continuous wavelets and frames on manifolds at the end of the introduction of 
[9], and for space considerations, will not repeat it here. In the last section of this article, we will present 
a detailed comparison of our methods with the methods of Narcowich, Petrushev and Ward [14], [15], 
who worked on the sphere; we will explain the similarities and differences between their methods and ours. 

To see how frames can be obtained from (7), suppose that, for any t > 0, K t is the Schwartz kernel of 
f(t 2 T). Thus, if F € L 2 (M), 

[f(t 2 T)F](x) = [ F(y)K t (x,y)d»(y), (8) 
Jm 

here /i is the measure on M arising from integration with respect to the volume form on M. Say now 
that J M F = 0, so that F = (I - P)F. By (7), 

A a (F,F) < (Y,\!\ 2 (a 2 ^T)FF) < B a (F, F). (9) 

3 

Thus 

A a (F,F) < J2(f(^T)F,f(a 2 ^T)F) < B a (F, F) , (10) 

3 

so that 

d»(x)<B a (F,F) (11) 

Now, pick b > 0, and for each j, write M as a disjoint union of measurable sets Ej^k with diameter at 
most ba? . Take xj.k € Ej.k- It is then reasonable to expect that, for any e > 0, if b is sufficiently small, 



Aa(FF) <Y,J J K ai (x,y)F(y)dfx(y) 
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and if Xj tk e Ej tk , then 



(A a - € )(F,F)<J2Y1 [ K °> y) F (y)Mv) 



2 



KE 1 ,k)<(B a + e)(F,F), (12) 



which means _ 

(A a e)(F,F) <J2J2\( F ^i-k)\ 2 < (Ba + e)(F,F), (13) 

j k 

where 

4>j,k{x) = fi(E^ k ) 1/2 K^(x^ k ,x), (14) 
or, formally, from (8), (since J(a 2 ^T) = f(a 23 T)*), 

tt> jtk {x) = KE 3 , k ) 1/2 [f(a 2 >T)5 X]k }(x). 

We will show that (13) indeed holds, provided the E^ k are also "not too small" (precisely, if they satisfy 
(16) directly below). In fact, in Theorem 2.4, we shall show (a more general form of) the following result: 

Theorem 1.1. Fix a > 1, and say c ,8 > 0. Suppose f e <S(M + ), and /(0) = 0. Suppose that the 
Daubechies condition (3) holds. Then there exists a constant Co > (depending only on M, /, a,c and 
So) as follows: 

For t > 0, let K t be the kernel of f(t 2 A). Say < b < 1. Suppose that, for each j e Z, we can write M 
as a finite disjoint union of measurable sets {Ej_ k : 1 < k < Nj}, where: 

the diameter of each E^ k is less than or equal to ba? , (15) 

and where: 

for each j with ba? < S 0} n(Ej tk ) > c (&a J )™. (16) 

(Such Ej yk exist provided cq and Sq are sufficiently small, independent of the values of a and b.) 

For 1 < k < Nj, define (f>j tk by (14-). Then if P denotes the projection in L 2 (M.) onto the space of 

constants, we have 

(A a -C b)(F,F) <^^|(f,^, fe )| 2 < (B a + C b)(F,F), 

j k 

for all F G (I - P)L 2 (M). In particular, if A a -C Q b > 0, then {4> jyk \ is a frame for (I - P) L 2 (M) , with 
frame bounds A a — Cob and B a + Cob. 

Thus, in these circumstances, if b is sufficiently small, {<j)j,k} is a frame, in fact a nearly tight frame, since 

f^~^ = l + 0(|(a-l) 2 (log| a -l|)|). 

Most of the work in this article will be in justifying the passage from (11) to (12). We will shortly sketch 
how this is done. First, however, we should explain in what sense the <f>j^ k look like a wavelet frame. We 
recall that, on the real line, wavelets are obtained through dyadic translations and dilations from one 
fixed function. 

Let us first, then, explain in what sense the passage from one j to another is analogous to dilation. For 
this, we need to understand very precisely the behavior of K t (x,y) near the diagonal. 
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If we were on R n , K t {x, y), the kernel of f{t 2 A), would be of the form t n ip{{x — y)/t) for some ip £ <S. 
(Here ijj = G, where G(£) = /(£ 2 )-) For any N, a, (3, there would thus exist CN,a,p such that 



t n+|a| + |;3| 



y 



N 



< c 



N,a,(3 



for all t, x, y. 

In Lemma 4.1 of [9], we showed that one has similar estimates on M, so that, in effect, K t {x, y) "behaves 
like" t~ n ip((x — y)/t) near the diagonal D (for some ip). Precisely, we showed: 

Lemma 1.2. For every pair of C°° differential operators X {in x) and Y {in y) on M, and for every 
integer N > 0, there exists Cn,x,y as follows. Suppose dcgX = j and deg Y = k. Then 



t 



n-\-j-\-k 



d{x,y) 



N 



XYK t {x,y) 



< Cn,x,y 



(17) 



for all t > and all x, y £ M. 



(Here d is geodesic distance on M.) In this sense, then, by (14), passage from one j to another (for the 
(j>j^k) is analogous to dilation. 

It remains to discuss if there is any sense in which the passage from one k to another (for fixed j) is 
analogous to translation. A general manifold M has nothing akin to translations, but in section 6 of [9], 
we discussed the situation in which M has a transitive group G of smooth metric isometries, and in that 
case there is a satisfying result. (Manifolds with such a group G are usually called homogeneous. Obvious 
examples are the sphere and the torus.) In [9], we observed that K t {Tx,Ty) = K t (x,y) for all T £ G 
and x,y <E M. Using (14), let us write <f>j t k — ^{Ej^Y^^PjM where tpj^{x) — K a j{xj^,x). Then, for any 
fixed j, we see at once that one can obtain all of the tpj t k by applying elements of the group G to any 
one of them. For example, on the sphere, for any fixed j, all of the tpj ^ are rotates of each other. This is 
a natural analogue for M of the usual requirement on R n that all elements of the frame at a particular 
scale j be translates of each other. 



Besides its pleasant and suggestive appearance, Lemma 1.2 (as stated above) is actually crucial, in en- 
abling us to rigorously pass from (11) to (12). Let us explain why. We fix a finite set V of real C°° vector 
fields on M, whose elements span the tangent space at each point. We let Vq — V U { the identity map}. 
Let Qj°(M) equal the set of smooth <p on M with integral zero. For x e M and t > 0, we define 

< 1 whenever y eM, 0<iV<n + 2 and Ye -p j 

For example, if we define <fl.{y) — K t {x,y), then by Lemma 1.2, there exists C > such that is in 
CM x j for all x e M and all t > 0. (The space M. x ,t is a variant of a space of molecules, as defined 
earlier in [11] and [13].) 

Our main technical tool will be the following fact. For each j £ Z, we write M as a finite disjoint union of 
sets Ej t k of diameter at most a J , and we choose points Xj t k £ Ej t k for all j, k. We then choose functions 
tfij,k,ipj,k € M Xj k , a i for all j, k. We then claim that the "summation operator" S defined by 



M 



X,t 



ip £ C™{M) : t n+dcsY 



d{x,y) 
t 



N 



Ytp(y) 



j k 



(18) 



G 



is bounded on L 2 (M), and that \\S\\ < C where C is independent of our choices of £y,fe, <Pj,ki V'j.fe- A 
precise statement is given in Theorem 2.3. To prove this result, we use the T(l) theorem to show that 
S is a Caldcron-Zygmund operator. (If M were W l , and if the ipj^ and ipj^ were obtained from fixed 
functions by means of translations and dilations in the usual manner, then this result is well known (see 
e.g. [11]).)' 

This result shows that, at least, the sum in (12) (or equivalcntly (13)) is finite, since it equals (SF,F), 
provided, in the formula (18), we set 

<Pj,k(y) = <Pi,k{y) = K ai(Xj,k,y)- (19) 

Once this is known, here is the spirit of how one can pass from (11) to (12). One can note that, by 
the fundamental theorem of calculus, the error made when replacing the outer integral in (11) by the 
Riemann sum in (12) can be expressed in terms of the integrals of appropriate derivatives over short 
intervals. Theorem 2.3 can then be used again to bound the integrand in this expression for the error, 
and the error itself is then small (in fact is 0(b)) because the intervals of integration are short. (It will 
actually be convenient to use identities such as b Jq H(bs)ds = J Q H(s)ds to convert the integrals over 
short intervals to integrals over [0, 1]; then the desired b will actually multiply the integrals.) The precise 
argument will be given in Theorem 2.4. 

A major advantage of working with functions of A is that one can now perform a space- frequency analysis. 
Based upon how well-localized a function F € L 2 is in space and in frequency, we can describe which 
terms in the summation 

j k 

are so small that they can be neglected. (Here, ipj^ is again as in (19).) In the analogous situation on the 
real line ([2], Chapter 2), one may do this by means of a time-frequency analysis. The analogue here is to 
do a "space-frequency" analysis (it is hardly reasonable to think of M as parametrizing time!). However, 
in contrast to the analysis for R in [2], we use the spectral theorem. We shall first do a frequency analysis. 
(Assuming, roughly, that F is well-localized in frequency, we describe which j can be neglected). This will 
reduce us to a finite sum. Then we will do a spatial analysis. (Assuming, roughly, that F is well-localized 
in space, we describe which fc can be neglected for each of our now finite set of j). 

Specifically, let = A < Ai < . . . be the eigenvalues of A on M, and let {u m : < to < oo} be an 
orthonormal basis of eigenvectors with Au m = A m w m for each to. If F = c m u m E L 2 (M.), and L > 0, 
let P[o,l]F = J2\ <l c m u m- Also let J > 1 be an integer, and suppose that F vanishes to order I > 1 at 
0. In Theorem 2.5 we shall obtain the estimate 

\\SF- J2 £m(^(^,*)MI2< ( C ob+^ti+-%) \\F\\ 2 + 2B a \\(I-P [0M )F\\ 2 , (20) 

n t , \ Qj d / 

j=-M k x 7 

where the constant c' L depends only on L, /, and a, the constant C'j depends only on J, /, a and M, and 
where B a is an upper frame bound for the frame {fi(Ej ^^ipj ^}. (c' L and C'j are identified in Theorem 
2.5.) 

Thus, say one wants to compute SF to a certain precision. One calculates the finite sum 

JV 

X! X^M^K^^fc)^ ( 21 ) 

j=-M k 
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for M, N sufficiently large; how large must one take them to be? One first chooses b, J, L, N so that the 
first, third and fourth terms on the right side of (20) are very small. Then one chooses M to make the 
second term on the right side of (20) very small as well. 

Let us explain, briefly, why one ought to expect a result like (20). Note that (13) says that, if F = (I—P)F, 
then 

(A a - e) (F, F) < (SF, F) < (B a + e) (F, F) . 
Thus, the arguments leading from (9) to (13) show that 



SF~ J2 mV 3 'A)F. 



In turn, by (7) and (4), 
that the truncated sum 



\,f\ (a? 3 A)F <~ CqF, where C = c/(2|logo|). Similar reasoning shows 



N N 
j=-M fe j=-M 



Let us then try to understand why, if F is well-localized in frequency, then -__ M |/| 2 (a 2j 'A)F <~ C F 
for M, N sufficiently large. That would then strongly suggest that we can truncate the series for SF with 
only a small error, as in (20). For A > 0, let 



5(A)- £ \f(a 2j \)f 



(22) 



the Daubechies sum, which is always between A a and B a , for any A, although (since / vanishes at and 
oo) the set of terms in the series which are sizable depends very much on A. Let 



N 



9m,nw= E i/(« 23 ' a )i 2 - 



-M 



If F - E m >i c m u m € L 2 (M), then 



(23) 



E m v j t)f -E c - 

j=-M 



m—1 



N 



E iffi^T)^ 

j=-M 



— E c m5M,Ar(A m )u ri 



m—1 



(24) 



If we had the Daubechies sum g here in place of gM,N, this would be a good approximation to CqF. In 
the actual situation, say F is well-localized, in the sense that all c m with X m > L are negligible. Then 
the rightmost member of (24) will be a good approximation to CqF provided 

ffM,Jv(A m ) — g{X m ) 

whenever Ai < A m < L. Recalling (22) and (23), and the fact that / vanishes at and vanishes rapidly at 
oo, we see that this happens provided that the interval [— M, N] contains all j with a? 3 X m ~ 1 whenever 
Ai < A m < L, or in other words, all j with C/L < a? 3 < C/X\. This will happen if M and N are 
sufficiently large. Thus, we expect that we can truncate with only a small error, as is precisely borne out 
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in (20). 

(20) reduces the approximate evalustion of SF to the calculation of the finite sum (21). It may well be 
possible to further reduce the number of terms that need to be considered, provided F is well-localized in 
space, e.g. if it is supported in a ball B of small radius. For, say Xj.k is far from this ball. By Lemma 1.2, 
for any / > 0, there exists Cj such that \fj,k(y)\ < Ci a ^~ n ^ d(xj t k, y)~ ! for all y e B, which is small; so 
\(F,ipj^)\ is small as well. In Theorem 2.6, we make these considerations precise, thereby performing a 
"spatial analysis" as well. 

Our frames should be viewed as a discretization of the continuous wavelets we presented in our earlier 
article [9]. However, it is important to notice that we will not actually make use of those continuous 
wavelets, since we get better results by discretizing Calderon's integral first (to obtain Daubechies' con- 
dition), and then discretizing as in (9) - (13). In this manner, we obtain frames if Daubechies' condition 
holds, and nearly tight frames as a — > 1 + , with nearly quadratic convergence of the ratio of the frame 
bounds to 1 - results that we would not obtain if we were to discretize our continuous wavelets in a more 
standard manner. (As an aside, we note that if one did want to discretize the consinuous 5-wavelets of 
[9] in a standard manner, one would want to require an additional condition such as (8) of that article, 
in order to be able to apply the fundamental theorem of calculus in the t variable. As we noted in [9], 
this additional condition does hold if K t is the Schwartz kernel of f{t 2 T).) 

In another article (already available [10]), we show that one can determine whether F is in a Besov space, 
solely from a knowledge of the size of its frame coefficients. In a future article, we hope to study the 
same question for Tricbcl-Lizorkin spaces. (The analogous problems on R" were solved in [6] and [7].) 

2 Frames 

We shall need the following basic facts, from Section 3 of [9], about M and its geodesic distance d. For 
x e M, we let B(x,r) denote the ball {y : d(x,y) < r}. 

Proposition 2.1. Cover M with a finite collection of open sets Ui (1 < i < I), such that the following 
properties hold for each i: 

(i) there exists a chart (Vj,0i) with Ui C Vi; and 

(ii) 4>i(Ui) is a ball in R n . 

Choose 6 > so that 35 is a Lebesgue number for the covering {Ui}. Then, there exist ci,c 2 > as 
follows: 

For any x € M, choose any Ui D B(x,3S). Then, in the coordinate system on Ui obtained from (pi, 

d(y,z)<c 2 \y-z\ (25) 

for all y,z € Ui; and 

ci\y-z\<d{y,z) (26) 

for all y, z £ B(x,S). 

We fix collections {Ui}, {Vi}, {4>i} and also <5 as in Proposition 2.1, once and for all. 
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• Notation as in Proposition 2.1, there exist c 3 ,c 4 > 0, such that, whenever x £ M and < r < 6, 

c 3 r n <^(B(x,r))<c 4 r n (27) 

• For any N > n there exists CV such that, for all x £ M and i > 0, 

/ [l + d(x,y)/t]- N dn(y)<C N t n (28) 
Jm 

• For any N > n there exists C' N such that, for all x £ M and i > 

■Y/ J.U — N 



f d(x,y)- N d»(y)<C' N t" 

Jd(x,y)>t 

For any N > n there exists C'^ such that for all x, y £ M and t > 0, 



(29) 



(30) 



M 

• For all M, t > 0, and for all E 1 C M with diameter less than Mt, if xo £ E, then one has that 
1 



M + 

for all x £ E and all y £ M. 



-[l + d(s,y)/t] < [l + d(a;o,i/)/t] < {M + !)[! + d{x,y)/t] 



(31) 



Now, fix a finite set V of real C°° vector fields on M, whose elements span the tangent space at each 
point. Using compactness, it is then easy to see that, if V is another such set, there is a constant C > 
such that for any / £ C 1 (M), and any x £ M, we have that 

max \X'f(x)\ < Cmax \Xf(x)\. 

As will soon be apparent, this implies that the choice of set V is immaterial (up to a constant) for the 
estimates which follow. For convenience, we choose V so that for every I with 1 < I < n, and every i, 
there is a vector field X in "P such that X = d/dxi on Ui (in the local coordinates obtained from (pi). 
(Notation as in Proposition 2.1). We also let 



Let 



Vo = V U { the identity map}. 
C5°(M) = jy> e C°°(M) : y ^ = oJ. 



For any x € M, we let 

M Xtt = LeC °°(M) : T+ dc s^ 



d(x, y) 



AT 



< 1 whenever y e M, 0<A^<n + 2 and Y £ V j 

(32) 
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Important Example Notation as in Lemma 1.2, for each x £ M, we define the functions 'p t x ,'^ t x on M 

by <^(y) = i^*(a;, y) and ip x (y) = K t (y, x). Then, by Lemma 1.2 there exists C > such that <p* and tp x 
arc in CM x j for all x £ M and all t > 0. 

Our main task in this section is to justify passing from (11) to (12). Our main technical tool will 
be the following fact. For each j £ Z, we write M as a finite disjoint union of sets Ej k of diameter at 
most a 3 , and we choose points Xj,k £ Ej^ for all j,k. We then choose functions <Pj,k,i>j,k G -M- Xj k , a i for 
all j, k. We then claim that the "summation operator" S defined by 

SF = E ^ E ^)(F, fiM,k (33) 

3 k 

is bounded on L 2 (M), and that ||5|| < C where C is independent of our choices of Ej t k, <Pj t k, 4>j,k- (A 

precise statement is given in Theorem 2.3 below.) 

If we recall the definition of the space M Xyt , in (32), we see: 

• For every C\ > there exists C 2 > such that whenever t > and d(x, y) < C\t, we have 

M v ,t C C 2 M x , t . (34) 

Say now x £ M, and R > 0. We say that a function w e C 1 (M) is a &itmp function for the ball 
B(x,R), provided supp u C B(x,R), \\cu\\oo < 1 and maxjegp HXwHoo < 1/i?. If w is a bump 
function for some ball B(x, R), we say that u> is an R-bump function. (This notion is modelled on 
the following situation on W 1 : take a C 1 function (p on R™ with support in the unit ball and with 
IMIc 1 < 1) an d look at o>(x) = (p(x/R).) Note also that even/ C 1 function on M is a multiple of 
an i?-bump function, where R — 2(diamM). 

• There is a constant C > such that for every R < 2(diamM), every i?-bump function lo, and every 
x, y £ M, one has 

\u(x)-w(y)\<Cd(x,y)/R. (35) 

To see this, we again use the notation of Proposition 2.1. Ifd(x,y) < 5, we choose Ui D B(x,3S), 
and use the mean value theorem in the local coordinates on Ui obtained from <j>i, to obtain (35). If, 
on the other hand, d(x, y) > 5, we simply note 

\lo(x) - u(y)\ < 2 - [2R/d(x, y)][d(x, y)/R] < [4(toanM)/S\[d(x, y)/R]. 

If T : C X {M) L 2 (M) is linear, we say that a linear operator T* : C^M) -> L 2 (M) is its formal adjoint 
if for all /, g £ C 1 we have 

<Tf,g>=<f,T*g>. 

T* is evidently unique if it exists. 

Recall that To — V U { the identity map}. We will be using the following form of the David-Journe T(l) 
theorem [4]. 

Theorem 2.2. There exist Co,N > 0, such that for any A > 0, we have the following. Whenever 
T : C£(M) -> L 2 (M) has a formal adjoint T* : C*(M) -» L 2 (M), and whenever T, T* satisfy: 

{i) \\Tw\\ 2 < AR n ' 2 and \\T*u\\ 2 < AR n / 2 for all R-bump functions uj; 

(ii) There is a kernel K(x,y), C 1 off the diagonal, such that if F £ C 1 , then for x outside the support 
ofF, (TF)(x) - / K(x,y)F(y)dy; and ' 
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(Hi) Whenever X (acting in the x variable) and Y (acting in the y variable) are in Vq, and at least one 
of X and Y is the identity map, we have 

\XYK(x,y)\ < Ad(x,y)- (n+dc z x+dceY *> 

for all x, y G M with x ^ y; and 

(iv) r(i) = r*(i) = o, 

then T extends to a bounded operator on L 2 (M), and \\T\\ < CqA. 

A simple proof of this theorem in the case M = W L may be found in [3] . That proof immediately adapts 
to our situation, in which M is a smooth compact oriented Riemannian manifold. 

Note: In verifying condition (i) of Theorem 2.2, one may assume that R < 2(diamM). For, set 
i?o = 2(diamM); it is evident that if R > Rq, then any i?-bump function cu is also an i?o-bump function. 
If we know that \\Tuj\\ 2 < AR n J 2 \ then we also know that ||Tw|| 2 < AR n l 2 . 

We can now prove the main facts about the summation operator (defined in (33)). 
Theorem 2.3. Fix a > 1. Then there exists C\,C 2 > as follows. 

For each j G Z, write M as a finite disjoint union of measurable subsets {Ej ^ : 1 < k < Nj}, each of 
diameter less than a? . For each j,k, select any Xj_ k £ Ej,k> an d select </?_,-. fe , Vj,fc with 
Vi,k,i>j,k € M-X-, k . a i- F° r F € C 1 (M), we claim that we may define 

SF = S Wjtk}t{ ^ k} F = ^2^2(i(E j!k )(F,ip jik )ip :jik . (36) 

j k 

(Here, and in similar equations below, the sum in k runs from k = 1 to k = Nj .) Indeed: 

(a) For any F G C 1 (M), the series defining SF converges absolutely, uniformly on M, In fact, if, for 
some R < 2(diamM), F = ui is an R-bump function, then the sum of the absolute values of the 
terms is less than or equal to C\ at each point o/M. 

(b) || SF || 2 < C 2 || F || 2 for all F G C\M). 

Consequently, S extends to be a bounded operator on L 2 (M), with norm less than or equal to C 2 - 
In fact, if we put T = S, then S satisfies the hypotheses of Theorem 2.2. 

(c) If F G L 2 (M), then 

SF = J2J2^ E ^)( F ^M,k 

j k 

where the series converges unconditionally. 

(d) If F,G G L 2 (M), then 

(SF, G) = J2Y1 (F, <Pj,k) W>j,k, G) , 

j k 

where the series converges absolutely. 
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Proof of (a). We first prove the second statement of (a). For each j, and each R < 2(diamM), we let 

C jtR = sup|(w,v?j,fc)| 

where the sup is taken over all i?-bump functions u) and all k. Since \{w,ipj tk )\ < ||w||i||<£j,fc||cx> we surely 
have that 



Cj, R < CR n a- ]n . 
On the other hand, we claim that if a- 5 < R, then 

Cj, R < Ca J /R. 

In fact, since J ipj^d/i = 0, we find 



\(u,Vj,k)\ 



< 



Lu(x)ip jik (x)d(j,{x) 
J (u(x) - u{xj >k )) <Pj,k{x)diJ,(x) 
J \uj(x) -u(x jlk ) || (pj. k (x)\d^(x) 

< C J R~ 1 d(x,x jtk ) \(fij. k {x)\dfi(x) 

< CV /i? 

by (28). Now, for any j and any y <G M, we have 

^2n{E jik ) \(u},tpj, k )\ \ipj,k(v)\ ^ C j,Rj2K E j,k) \ipj,k(y)\ 

k k 

< CC j>R a-* n ]T »(E jlk ) [1 + d(y, Xi , k )/a? 



-n-1 



< CC 3 . R a-> n / [l + d(y,x)/a j 



dfj,(x) 



again by (28). (In passing from (42) to (43), we used (31).) 
Taking the sum over j we have: 



j k 



< c 
<c 1: 



ai<R 



ai>R 



(37) 
(38) 

(39) 



(40) 

(41) 
(42) 

(43) 
(44) 



(45) 
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by (37) and (38). This proves the second statement in (a). Since every F <G C 1 is a multiple of an i?-bump 
function, with R = 2(diamM), the first statement in (a) is a consequence of the two calculations above, 
which ended with (44) and (45). □ 



Proof of (b). To begin, let us set 

K,t= $eC°°(M) : t n 



d(x,y) 



N 



<%) 



< 1 whenever y <G M, < N <n + 2 



Claim: There exists C 3 > (independent of our choice of the E^ k and Xj^) as follows: Suppose {Qj.k} 
and {*j,fc} are two systems of functions such that 3>j,fc, ^j,k & N x . k , a i for all j, k, and suppose < J < 1. 
Then 

Kj{x,y) : "V/^i \*j,k(x)*j,k(v)\ < C 3 d(x,y)- n - J 

j k 

for any ije M, x ^ y. 

To prove the claim, note that for any j, 

» iJ n( E i,k) \*j,k(*)*j,k(v)\ < a- i(J+2r,) EM£i, fe ) [1 + d(x,^, fc )/^r ( " +2) [1 + fc )/a']~ (n+2) 
fe k 

< Ca^ J+2n ) [ [1 + d(x, z) /a'] " ( " +2) [1 + d(y, z) /a'] ~ ( " +2) d/x(s) 

JM 



< Ca-^ J+ ") [l + d{x,y)/a j 



-(n+2) 



by (30). 

Say x ^ y , and let j be the integer satisfying o- 70 < d(x,y) < a JO+1 . Then we have 



d(a^) J+ "Mz,</) < C^a-i^dfi,!,)^!!!^,!/)^] 



ji-n-2 



< c 



< c* 



< c. 



3(J+n) n ]a(J+n) + \~~ „ jt.l m^j^.l n I „ ( „ +2 ) I., -.„, ) 



^ a -Oo-j)(J+n) + ^ a (jo-j)(J-2) 
J>jo 3<jo 



since we are assuming that < J < 1. This establishes the claim. 

We return now to (36). By definition of A^t, a function ip E Ai x ,t if and only if t dcsY Yip G Af x ,t 
for every 7 e Po- We apply this with t — a? (j G Z). It therefore follows from the "claim" that, 
whenever X (acting in the x variable) and Y (acting in the y variable) are in V , then 



£5>(^,fc) \X^, k (x)YJp jtk (y)\ < C 3 d(x iy )-^ +dc z x + d ^ Y \ 

j k 



(46) 
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This has a number of immediate consequences. Taking X and Y to be identity maps, we now see that 

j k 

is the kernel of 5, in the sense that if F e C 1 , then 



[SF](x) = f K{x,y)F(y)dn(y) 
Jm 



whenever x is outside the support of F. (This is because, for fixed x, the series defining K converges 
absolutely for y in the support of F, and the sum of the absolute value of the terms is bounded inde- 
pendently of y.) Secondly, since we could, in local coordinates on a Ui, take X in (46) to be any d/dxi 
(and Y to be any d/dyi), we see that K(x,y) is C 1 off the diagonal. Finally we see that we may bring 
derivatives past the summation sign in (47), and thus, for any X, Y e Vo, 

\XYK{x,y)\ < C 3 d(x,y)-( n+dcsX+dcsY ^ 

for x ^ y. 

This shows that 5 satisfies conditions (ii) and (Hi) of Theorem 2.2; we need to show it satisfies the other 
conditions. By (a), it is evident that if F, G £ C 1 (M), then 

(SF, G) = J2Y1 ^ E ^)(F, <Pj,k)(1>i,k, G), 
j fe 

where the sum converges absolutely. Consequently, the formal adjoint of S is S*, where, if G € C 1 , 

j k 

Surely 51 — 5*1 — 0, since J ipj^ — J ipj^ — for all j, k. 

That leaves only condition (i) of Theorem 2.2. Suppose then that uj is a bump function for the ball 
B(x ,R), where R < 2(diamM). If d(x,x ) > 2R, then 



\[Sw](x)\ < J \K(x,y)u(v)\dv(y) 

<C [ d(x,y)' n dy 

Jd(y,x )<R 

<CR n max d(x,y)~ n 

d(y,x )<R 

< CR n d(x,x )- n . 

Recalling (a), we now see that 

\\Su\\l= ( \[Scj](x)\ 2 dn(x) + f \[Su;}(x)\ 2 dKx) 

J d(x,x )<2R Jd(x,x )>2R 

<C (2R) n + R 2n j d(x,x )- 2n dti(x) 

J d(x,x )>2R. 

< C[R n + R 2n R- n ] = CR n , 
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by (29). Hence 

\\Su\\ 2 < CR n/2 , 

This establishes condition (i) of Theorem 2.2, and completes the proof of part (b). □ 

Proof of (c). It is evident, by (a), that (c) holds for F e C^M). 
Suppose now that T C 1? is finite, and define : L 2 — > C 1 by 

^= I] l*(Ej,k){F,<Pj,k)il>j,k 

By (b), 5^ : L 2 — > L 2 is bounded, with norm HS^H < C 2 for all N. (Indeed, in the formula for S, we 
have just replaced the ipj. k by the zero function if j, k £ T , and the zero function is surely in all Ai x . k a j .) 
Since (c) is true for F £ C 1 , which is dense in L 2 , it now follows for all F e L 2 . □ 

Proof of (d). This follows at once from (c), since a series of complex numbers conveges absolutely if and 
only if it converges unconditionally. 

This completes the proof of Theorem 2.3. □ 

We are now almost ready to justify the transition from (11) to (12). In order to do so, we need to choose 
the Ej,k rather carefully. We will need to utilize the following simple fact: 

• Say t > 0. Then there exists a finite covering of M by disjoint measurable sets E\, . . . , En, such 
that whenever 1 < k < N, there is a y k € M with B(y k ,t) C E k C B(y k ,2t). Thus, by (27), there 
is a constant c' , depending only on M, such that ^i(Ek) > c' Q (2t) n , provided t < 8. 

To see this, choose a maximal collection of disjoint balls {B{y kl t) : 1 < k < N}, set B k = B(y k , t), 
and set B' k = B(y k , 2t). Since the collection is maximal, each point y G M is at distance less than 
t from some point in U k=1 B k , so the distance from y to some y k is less than 2t. Therefore the 
B' k cover M, and we need only set E\ = B[ \ (B 2 U . . . U Bjy), and recursively, for 2 < k < N, 
Ek = B' k \ {Ei U . . . E k -i U -Bfe+i U . . . B N ). 

We conclude: 

• Say c' Q is as in the last bullet point, and say < b < 1. Then, for each j € Z with 6a- 7 < 25, we can 
write M as a finite disjoint union of measurable sets {Ej k }, where the diameter of each Ej tk is less 
than or equal to ba? , and where, for each j with 6a- 7 < 28, n(Ej^ k ) > Co(ba : >) n . 

This follows from the last bullet point, with t = ba^ /2. 

The following theorem now justifies the transition from (11) to (12), and in particular implies Theorem 
1.1. 

Theorem 2.4. (a) Fix a > 1, and say co,Sq > 0. Suppose f e <S>(M + ), and /(0) = 0. Then there exists 
a constant Co > (depending only on M, /, a,c and So) as follows: 

Let J C Z be such that for some N , j <G J whenever \ j\ > N. (We will often take J = Z.) For t > 0, 
let K t be the kernel of f(t 2 A). For 1,1/6 M, set 

®x,t(y) =K t (x,y). 



1G 



Say < b < 1. Suppose that, for each j e J , we can write M as a finite disjoint union of measurable 
sets {Ej^k ■ 1 < k < Nj}, where: 

the diameter of each Ej^ is less than or equal to ba? , (48) 

and where: 

for each j with ba? < So, ii(Ej^) > Co(ba 3 ) n . (49) 

(By the last bullet point, we can surely do this if c < c and 5q < 25.) Select Xj^ G Ej ^ for each 
j,k. By Lemma 1.2, there is a constant C [independent of the choice of b or the Ej t k), such that 
. G CM Xj k , a j for all j,k. Thus, if for 1 < k < Nj, we set 



. ., if j G J , 
otherwise, 



we may thus form the summation operator S J with 

j k 

and Theorem 2.3 applies. 

Let Q J = J2jeJ \f\ 2 ( a2:iT ) (strong limit, as guaranteed by (7).) Then for all F e L 2 (M), 

\((Q J -S J )F,F)\<C b(F,F) (50) 
(or, equivalently, since Q J — S J is self-adjoint, \\Q J — S J \\ < Cob.) 

(b) In (a), take J — Z; set Q = Q 1 , S ~ S z . Use the notation of (7), where now T = A on L 2 (M.); 
suppose in particular that the Daubechies condition (3) holds. Then, if P denotes the projection in L 2 (M.) 
onto the space of constants, we have 

(A a - C b)(I -P)<S<(B a + C b)(I - P) (51) 

as operators on L 2 (M). Thus, for any F E (I — P)L 2 (M), 

(A a - Cob)\\F\\ 2 <J2KE 3 ,k)\(F, ^,k)\ 2 < (B a + Cob)\\F\\ 2 , 

so that, if A a — Cob > 0, then {^(Ej.kY^^Pj.k} • k is a frame for (I — P)L 2 (M.), with frame bounds 
A a — Cob and B a + C b. 

Note By (4), B a /A a = 1 + 0(\(a - l) 2 (log \a - 1|)|); evidently (B a + C a b)/(A a - Cob) can be made 
arbitrarily close to B a /A a by choosing b sufficiently small. So we have constructed "nearly tight" frames 
for (I -P)L 2 (M). 

Proof. We prove (a) . To simplify the notation, in the proof of (a) , j will always implicitly be restricted 
to lie in J, unless otherwise explicitly stated. 

Since So only occurs in (49), we may assume So < S; for otherwise, we may replace So by S, and (49) still 
holds. 
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Since Q J and S J are bounded operators, we need only show (50) for the dense subspace C 1 (M). For 
b > 0, we let = log Q ((5 /&), so that j < is equivalent to 6o J < <5o- Observe that 



\((Q J ~S J )F,F)\ 



where 



and 



j k 



<I+II, 



E 



j>n b k 



For 77, we need only recall that, for any L > 0, 

lim = in C°°(M x M). 



(52) 



(This follows at once from the eigenfunction expansion of K t ; see Section 4 of [9], especially (25) of that 
article, and the comments that directly follow it.) In particular, there exists C > such that, for any x 
and any t > 1, ||3> x ,t||2 ^ CA- Accordingly, 



II < C E a->(M)||F||2 < Ca- n iF||2 = 

j>f2 b 



Thus we may focus our attention on I. We have 



EE / K^.^fw-EM^)!^,.^)! 

E E / [K^,^)| 2 - |(^^,F)| 2 ] d»(x) 
j<n b k Je o," 



(53) 



For each (j,k) with j < fib and 1 < k < Nj, we select ij.k such that B(xj.k,3S) C l/^ fc . (Since then 
diamiJ^fc < feci-' < <5, we then also have that .E^fc C B(xj t k,S) C [7^ fc .) For each i with 1 < i < AT, we 
let Si denote the set of (j, k), with 1 < k < Nj, such that ij^ = i- The summation in (53) may then be 
written as J2i=i ^ where 



E 



(j,k)eSi " a i- k 



{®x,ai,F)\ -\{<&x hk ,ai,F)\ dfl(x). 



(54) 



For (a), we only need show that each |Jj| < C 1 1 1 1 2 - 

Let us then fix i and use local coordinates on Ui. In these coordinates, Ui is a ball in W 1 . On Ui, we 
may write dji{x) = h(x)dx, where h is smooth, positive and bounded above. By (26), if (j, k) £ Si, then 
Ej,k Q B(xj^,ba : >) C jy e 1" : \x — Xj t k\ < ba^ /c\). In the integrand in (54), we may therefore write 
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x = x^ k + (ba 3 /ci)w for some w in the open unit ball in R n , which we denote by B. Changing variables 
in (54), we now see that 



F )\ \ [XE ilh h] (x)dx 



(ba 3 /cj) n [\(®x j , k +(bai/c 1 )w, a i,F)\ 



\{^ Xj , k ,ai,F)\ H jtk {w)dw, 



(55) 



where 



H j,k( w ) = iXE Jtk h](x^ k + (ba 3 /a)w). 

(This is interpreted as zero if Xj, k + (ba 3 /c\)w ^ Ej ik . Moreover the integrand in (55) is interpreted as 
zero at any point w where Hj yk (w) — 0.) Thus 

h= v{Ej,k) I \(^ Xj , k +(boJ/c 1 )w,ai,F)\ - \(® Xjtktai ,F)\ G jjk (w)dw, (56) 



where 



(ba 3 ) 7 



-H jik {w) 



Ci/j,(Ej tk ) 3 ' 1 

Note that, by (49), there is a constant C such that < Gj yk (w) < C for all (j,k) € 5, and all w € B. 

(Again, the integrand in (56) is interpreted as zero at any point w where Gj <k (w) — 0.) 

Applying the fundamental theorem of calculus, and recalling that Ui (being a ball) is convex, we now see 

that 



| (^Xj, fc +(6aJ / Cl )sw,ai > F}\ ds G Jyk (w)dw. 



(57) 



U,k)eSi 

For (j, k) e w e B and < s < 1, y e M, let us set 



J,k(y) = \J G i,k{w)<$> x . k+(ba i/ c ^ )sw , a i (y) 

= \jGj,k(w)K aj (x jik + (ba 3 /a)sw,y) 
(interpreted at zero if Gj. k (w) — 0.) Also note that, for y e M, 



/ (j I o 

\JG jik {w) — $ x . k+{baj/ci)sWtai {y) = y/Gj, k {w)—K aj (x jik + (ba 3 / Cl )sw,y) 



= Hj; k s (y) 



where 



i>7; k s (y) 



-i n 



m— 1 



,dK„ 



dx r 



(xj,k + (ba 3 /a)sw,y) 



Now, if w is such that Gj. k (w) =^ 0, then x^ k + (ba 3 /c\)w G Ej, k C Ui, so Xj ik + (ba 3 /c\)sw G Ui for all 
< s < 1. Since \[xj yk + (ba 3 /ci)sw] — xj ik \ < a? jc\, by (25), we have that d(xj yk + (ba 3 /ci)sw,Xj tk ) < 
(c2/ci)a 3 . By (17) and (34), there is consequently a C > such that 
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for all (j, k) e Si, w e B and < s < 1. Now let us set ip J' k s , ipj k s = whenever (j, k) £ Si (here j runs 
from — oo to oo, and k runs from 1 to Nj). We now find from (57) that 



Ii = b KEj,k) f f 



(j,k)ES. 



dsdw 



dsdw 



= b IJ! [ <s w>w }F > F) + ^OTW^ 

<Cb\\Ff 



by Theorem 2.3. (The interchange of order of summation and integration is justified by the dominated 

convergence theorem and the second sentence of Theorem 2.3 (a).) This proves (a). 

To prove (b), we need only show (51). But from (50) and (7), we have that, if F = (I - P)F e L 2 (M), 

then 

(A a - C a b)\\F\\ 2 < (QF,F)-C b\\F\\ 2 < (SF,F) < (QF, F) + C b\\F\\ 2 < (B a + C b)\\F\\ 2 . 

If, on the other hand, F e L 2 (M) is general, we have SF = S(I — P)F, since all (pj,k have integral zero. 
Since 5 is self-adjoint, (SF, F) = (S(I - P)F, (I - P)F), so in general 



(A a - C b)\\ (I - P)F\\ 2 < (SF, F) < (B a + C b)\\ (I - P)F\f 



as desired. 



□ 



In Theorem 2.4 (b) {/j,(Ej^ k ) 1 / 2 tpj_ k } is a frame. This is, of course, an infinite set, and so for practical 
purposes, we must explain which terms in the summation 



SF = v( E j,k)(F, <Pj,k)<Pi,k 



(58) 



are so small that they can be neglected. As we stated in the introduction, our plan is first to do a 
frequency analysis. (Assuming, roughly, that F is well-localized in frequency, we describe which j can be 
neglected). This will reduce us to a finite sum. Then we will do a spatial analysis. (Assuming, roughly, 
that F is well-localized in space, we describe which k can be neglected for each of our now finite set of 
j)- 



2.1 Frequency Analysis 

We begin with some motivation for the result we seek. 

Even without using frequency analysis, it is not difficult to see that the terms with j large and positive 
will contribute little to the sum (58). Indeed, for any L > 0, there exists cl > such that if F e L 2 , then 

Cj.L := \(F,tp jik )\ < ||F|| 2 ||^ ;fe || 2 < c L a- Lj \\F\\ 2 , 

by (52). Following (41) through (44) with Cj t L in place of C^r, we see that for each j, 

5>(^,fc) \(F,Vj,k)\ \<PjAv)\ < CC 3+- (59) 
k 

Taking the sum over j > N, say, we have: 

]T 5>(25i,fc) \(F, <p jtk )\ \^,k(y)\ <CJ2 C jlL < C L \\F\\ 2 a- LN , (60) 

j>N k j>N 
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which goes to zero rapidly as N — > oo. Here we need no information about F except its L 2 norm. 

On the other hand, it is not hard to see that if F has a little smoothness, then the terms with —j 
large and positive will contribute little to the sum (58). If, for instance, F is Holder continuous of ex- 
ponent a G (0,1), the arguments of (39) through (40) show that \{F,ipj t k)\ < Ca ja , so if we argue as in 
(59) and (60) above, we see that 

]T 5>(^,fc) \(F,<Pj,k)\ \<Pj,k(y)\<C E a ja <Ca~ Ma , 

j<-M k j<-M 

which also goes to zero rapidly as M — > oo. 

We leave the detailed study (through frames) of smoothness spaces to a later article ([10]). Here we 
instead look at the terms, with \j\ large, through frequency analysis. 

Specifically, say F — a i u i is the expansion of F with respect to the orthonormal basis of eigenfunctions 
{ui}. We expect that some knowledge of the decay of the {a;} should imply something about the decay 
of the terms in (58) as j —> — oo. After all, for instance, if the a; decay quickly enough, then F is C . 
Our next result gives us such information, in considerable generality. Moreover, the fact that the terms 
in (58) become negligible as j — > oo, will be seen to be a consequence of the fact that, if to > 1, then the 
eigenvalues of A satisfy \ m > Ai > 0. The proof of the result will depend on Lemma 2.3 of our earlier 
article [9]. 

Theorem 2.5. In the situation of Theorem 2.4 (b), suppose in fact that f has the form f(s) — s fo(s) for 
some integer I > 1 and f € <S(M + ), (/o ^ 0). Suppose J > 1 is an integer, and let Mj = max r>0 \r J f(r)\. 
For any L > 0, we consider the spectral projectors P[q,l] (so that, if F = ^a m u m e L 2 , then P[q.l]F = 
J2x m <L a "i u m)- Then if F e L 2 , and any M, N > 0, we have 

N d C 

\\SF- Y, EM^i,fc)<^yj,fc>Vi,fcll2 < (C 6+ + -4^)11^112 + 2 J B„||(/-Pto,£])i ; '||2, (61) 

H I 1 il il 

j = -M k 

where c' L = (L 21 ||/ |lL)/(« 4 ' - 1), and C = Mj/[(a iJ - l)Af ]. 



Proof. First note, that in proving this result, we may assume F = (I — P)F; otherwise, replace F by 
(I — P)F (this will not affect the left side of (61), and will not increase the right side). It then follows 
that (I-P [Q , L] )F= (I-P [XlM )F. 

Let J = {j e Z : j < —M or j > N}. Then the left side of (61) is simply \\S J F\\ 2 , in the notation of 
Theorem 2.4 (a). By that result, the left side of (61) is less than or equal to Co&||F|| 2 + HQ^-^lb- But in 
the notation of Lemma 2.3 of [9], Q J — h(A) — /im,jv(A). Thus the theorem now follows at once from 
Lemma 2.3 of [9] (with r\ in that lemma equalling Ai). This completes the proof. □ 

Thus, say one wants to compute SF to a certain precision. One calculates the finite sum 

JV 

Y ^Z^ E 3^)(F^ ] ,k)^P 3 .k (62) 

j=-M k 

for M, N sufficiently large; how large must one take them to be? One first chooses b, J, L, N so that the 
first, third and fourth terms on the right side of (61) arc very small. Then one chooses M to make the 
second term on the right side of (61) very small as well. 
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We emphasize that C'j in (61) does not depend on L, so that the negligibility of the terms as j — > oo 
depends only on ||-F||2, as we observed above. C of course depends on /, a and also on the manifold M 
(more specifically, on the least positive eigenvalue Ai). 



2.2 Spatial Analysis 

The summation in (62) is a finite sum, but, as we have motivated in the Introduction, if we only want 
to know this summation to within a certain precision, then it may well be possible to further reduce the 
number of terms that need to be considered, provided F is well- localized in space. 

In order to carry out the analysis, it is convenient to first make a general remark. Let 1 C {(j, fc) : 1 < 
k < Nj} be arbitrary. Let 

SxF= Y »( E 3,k){F,yj,k)Vj,k- (63) 
G\fc)ex 

Then Si is a positive operator on L 2 (M.); let \/Sx be its positive square root. For any F e L 2 , 

\\V&F\\l = (S X F,F) = Y, KE 3 ,k)\(F,^,k}\ 2 < B\\F\\l 
i, feel 

where B is an upper frame bound for the frame {^{Ej^Y^tPj.k}- Accordingly, ||\/5r|| < V^B, so 

\\Sx\\<B. (64) 

Also, for any F e L 2 , \\S X F\\ 2 < \\y/Si\\\\^F\\ 2 , so 

\\S X F\\ 2 2 <B(S X F,F). (65) 

Thus S X F will have small L 2 norm provided (SjF, F) is sufficiently small. Accordingly, the terms in (62) 
corresponding to (j, fc) e 1 can be neglected, provided (SxF, F) is sufficiently small. 

We may now prove out result on spatial analysis: 

Theorem 2.6. In the situation of Theorem 2.4 (b), say I > 0. Then there exists a constant C\ > 
(depending only on I,f and M), as follows. 

Let B be an upper frame bound for the frame ^{Ej^Y^^Pj.k}- Say M, N > 0, x <E M, and R > 0. Let 
X be the characteristic function of a setT C M. For each j with — M < j < N, let Cj > be a constant, 
and let 

1 = {(j, fc) : —M < j < N, 1 < k < N v d{x oM ,T) < ( Cj + l)a>} 

Then 

N 

Y X]M(-Ej,fc)(^^j,fe)^,fe - Y ^ E 3,k){ F ^j,k)fi,k 
'=— M k j,k£2 

1/2 



<C\B^ 2 



N 



M (L) Y a ^" c " 
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-M 



\\ X F\\ 2 + B\\(l-x)F\\ 2 



(66) 



for all F e L 2 (M). 
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Proof. Let T = {(j,k) : -M < j < N, 1 < k < Nj, d(x jik ,T) > ( Cj + l)a j }; then, in the notation of 
(63), the left side of (66) is \\Si>F\\ 2 . Surely 

\\Sx>F\\ 2 < \\Sx>(xF)h + \\Sx>([l - x}F)h < IIMx^lh + B||(l - x)F\\ 2 (67) 
by (64). On the other hand, if (j, k) € I' , then 



\(xF,V>j,k)f < 



\( X F)(y)\\^Ay)\My) 



< 



\<Pj,k(y)\ Mv) 



WxF\\ 



(l + d(x j!k ,y)/a J ) 21 dji{y) 



For each j with — M < j < N, let Ej = U -E-j.fc- If £ € Ej, then x G Ej, k for some (j, k), so for 

{fc: (j,k)el'} 

each y G T, d(x,y) > d(xj. k ,y) — d(xj ik ,x) > Cja^ . Thus, by (31) and then (29), we see that 



(Sx'(xF)AxF)) = £ M^vOKx^,*)! 5 

(j,fc)ex' 



< 



c? / E a ~ 2in i i (E j ,k) (i + fe, y)/« J ) 27 d/i(») llx^lll 



(j',fe)ex 

N 



j=-M 
N 



< C 2 2 21 ]T a- 2 -" 1 / / (l + rf^y)/^) "^(^(j,) 



r je 



< C 2 2 2/ ^ a- 2 -*" / / {d(x,y)/a*) 21 d^(x)d^y) \\ X F\\ 2 2 

i—-M J {x: d(x,y)>Cjai} 



j=-M 
N 



< C'jCp 21 £ «" J "c;- 2/ / d/x(y) ||xF|| 

j=-M "' r 
AT 

< c;c 2 2 2/ ^ a-^ Cj "- 2I M (r) Hx^lll 

j=-M 



Accordingly, by (65) we find 

WSx^h^C^B^T)) 1 ' 2 



N 
j=-M 



1-21 



1/2 



\\XF\\ 



2- 



with C\ depending only on /, / and M. If we combine this with (67), we see that the proof is complete. □ 

In all, say one wants to compute SF to a certain precision. By using frequency analysis, one can reduce 
the problem to computing X^-mS/c l I (Fj tk ){F,(pj, k ) L Pj,k to a certain precision. If F is well-localized, 
by Theorem 2.6, one can further reduce the problem to computing a sum J^. keX fx(Ej jk )(F, <Pj,k)fj,k, 
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by first choosing the set T so that B\\(I — x)_F|| 2 is small, and then choosing the numbers Cj so that 
SjL-M a _:, "c™ _2/ is sufficiently small. This will entail making the numbers Cj sufficiently large. The 
needed computations will increase as j gets smaller, since each a _jn c™ _2/ needs to be small (forcing 
larger requirements on cj as j gets smaller) , and since the number of k with (j, k) € X also increases as j 
gets smaller. Thus it is best to use frequency analysis to make M as small as possible. This will be most 
practical if F has a little smoothness. 



3 Needlets and Mexican Needlets on the Sphere 

In order to implement our frames, one of course would need to calculate approximately the frame elements 

4>j,k = K E j,k) 1/2 K aj (x jtk ,x), 

for some collection of Ej^ satisfying (48), (49) (for any given b e (0, 1)) and for some choice of xj.k € Ej^. 
In Section 6 of [9], we focused on the cases where M is the torus T 2 or the sphere S 2 . In these cases it is 
not difficult to explicitly choose suitable Ej_k- Moreover, in Section 6 of [9], we gave explicit approximate 
formulas for K t (x,y), if f(s) = se~ s (the Mexican hat case), for the torus and for the sphere, which 
therefore should make implementation feasible. In each case (T 2 or S 2 ), we gave two approximate 
formulas for K t (x, y), one which converges quickly for large t, and one which converges quickly for small 
t. In this section we focus on the sphere S 2 , where related work has been done by Narcowich, Petrushev 
and Ward ([14], [15]). We shall compare our methods with theirs, in detail. 

Specifically, on S 2 , in [9] we noted that K t (Tx,Ty) — K t (x,y) for any orthogonal transformation T. 
From this we concluded that K t (x,y) was actually a function of x ■ y, say 

h t (x ■ y) = K t (x,y). 

From the eigenfunction expansion of K t , we found that, for —1 < y\ < 1, if y = (t/i, . . . , ys) £ S 2 , 

oo 

4nh t ( yi ) = K t (N,y) = ^(2Z + l)f(t 2 l(l + l))P ; 1/2 (yi). (68) 
z=o 

1/2 

Here P l is a Gegenbauer polynomial. This series converges quickly when t is large. When t is small, 
and f(s) = se~ s , we obtained the approximation 

-e 2 /it 2 6 2 

47rh t (cos0) p — [(1 - ^)p(t,0) -t 2 q(t,6)], (69) 

where 

, m „ t 2 t 4 At & t 8 9 2 A 2t 2 4t 4 At\ 
p(M) = 1+ 3 + T5 + 3T5 + 315 + 4 ( 3 + 15 + 105 + 315 } 

and 

, ^ 1 2t 2 At A At 6 6 2 , 2 8t 2 At 4 , 
9(t ' fl) = 3 + 15 + i05 + 315 + 4 ( 15 + lO5 + W5 ) 
Maple says that when t = .1, the error in the approximation (69) is never more than 9.5 x 1CP 4 for 
any 9 e [— ir, n], even though both sides have a maximum of about 100. (To obtain rigorous bounds on 
the error is research in progress, which we expect to complete soon.) One can derive similar types of 
approximations if instead f(s) — s r e~ s for any integer r > 1, by applying repeated t derivatives to the 
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above formulas. (See the remark at the end of Section 6 of [9].) If in (69) we approximate p <~ 1 and 
q <~ 0, we would obtain the formula for the usual Mexican hat wavelet on the real line, as a function of 
9. This is to be expected, since on W 1 , the Mexican hat wavelet is a multiple of Ae~ A 5, the Laplacian 
of a Gaussian, the function whose Fourier transform is |£| 2 e~^' . Figure 1 is a graph, obtained by using 
Maple, of 4-7T h t (cos6), for t = 0.1, with 8 going from — n to tt on the horizontal axis. 
The related frames of Narcowich, Petrushev and Ward have been dubbed "needlets" , and used by astro- 
physicists to study cosmic microwave background radiation. (See, for instance, [1], [12] and the references 
therein.) We will therefore call the frame we obtain when f(s) — se~ s (or, more generally, f(s) = s r e~ s ), 
by the name Mexican needlets. The formulas above suggest that Mexican needlets have strong Gaussian 
decay at each scale, and that (at least for small r) they do not oscillate to an extent that would make 
implementation, directly on the sphere, prohibitive. 

We now give a detailed comparison of needlets and Mexican needlets; there are advantages and disad- 
vantages in each approach. In this discussion, H will denote L 2 (S n ) and Pi will denote the projection 
onto Hi, the space of spherical harmonics of degree I. 

In [14], [15], the authors considered only M = S n . In place of our / they consider only smooth g with 
compact support within [1/2,2] (we shall henceforth call such an g a "cutoff function"). They also ob- 
tained their frames from the kernel K t of an operator, K t (N, y) being given by a formula similar to our 
(68), but with g(tl) in place of f(t 2 l(l + n — 1)). (Thus they were not actually considering functions 
of t 2 A, or equivalently functions of t\/A, but rather functions of tA4, where M. is the first-order pseu- 
dodifferential operator E; ZP;. This is a minor distinction, however.) As we shall explain, the principal 
advantages of using cutoff functions is that the authors are then able to obtain tight frames, and that the 
frame elements on non-adjacent scales are orthogonal. The principal disadvantages, as we shall explain, 
is that there is no reason to expect explicit formulas on the sphere, Gaussian decay at each scale, or lack 
of oscillation for needlets. 

Let us begin our detailed discussion of needlets by explaining how Narcowich, Petrushev and Ward obtain 
a tight frame by using cutoff functions: 

1. One can easily choose a cutoff function g so that Ejl-oo lff(2 J " s )| 2 = 1 f° r all s > instead of this 
being only approximately true (recall our (3) and (4)). In fact, if g is a cutoff function, there are 
only two nonvanishing terms in the sum for any given s > 0. 

2. For any m > 0, let V m — YlT^i Pi- One can choose a finite set of cubature points {x m j} and 
positive numbers X m ,i such that for any F £ P m Ti, one has that J s „ F = Ej X m ,iF(x m ^). Thus 
when evaluating the integral of a function known to be in V m H, one can use cubature to evaluate 
the integral instead of approximating the integral by a Riemann sum (as we did in passing from 
(11) to (12)). 

One can construct tight frames out of such g, the plan being as follows. Suppose {fj}^ is a sequence 
of positive real numbers which increases to oo, with no worse than polynomial growth. Let M. be the 
self-adjoint operator Ez*li v iPi- ( in [14], [15], the authors take M. = E^ as in the last paragraph. We 
instead would look at M = y/A = J2 \AC + n— l)Pj.) For t > 0, let K t (x,y) be the kernel of g(tM). 
(We would choose g(s) = /(s 2 ), so that g{ty/A) = f{t 2 A).) Since E^L-oo lf( 2Js )l 2 = h for all s > 0, we 
have E^L-oo \g( 2i M)\ 2 = I - P . Thus, if F e {I - P Q )H, we find (analogously to (9), (10)) that 



OO OO 

F|| 2 = <£ | 5 (2^M)| 2 F,F> = E h{*M)F\\l 



(70) 
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Figure 1: 47r/i t (cos 6) on S 2 for t = 0.1 
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But since g is supported in [1/2,2], we may, for each j, choose with g(2^vi) — whenever I > so 
that g(2 j M) : TL — ► Vuj^H. As is well known, for any Z, the product of two elements of P/H is in ViiH. 
Accordingly \g(2^M)F\ 2 e V 2 i{j)H, and so || 3 (2JA1)F||| = / |.g(2-?.M)F| 2 may be evaluated exactly by 
cubature (instead of our having to approximate it by a Riemann sum, as in (12)). We find that 

00 00 00 

\\f\\1= J2 \\g(2 j M)F\\t= J2 Y,^mM^M)F\\x 2l{jU )= J2 £K F A*>I 2 < 

j = — OO j — — OO « J — — OO 2 

where 

<f>i,i{v) = \J X 2i{j),i K2i(x jti ,y) (71 ) 

(analogously to (14)). The {4>j,i} are therefore a normalized tight frame, not just a nearly tight frame. 
Moreover, the constraints on the supports of the dyadic dilates of g, easily imply that frame elements at 
non-adjacent scales are orthogonal. 

We turn to apparent disadvantages of needlcts. In the study of CMB, they are not implemented directly 
on the sphere, for several reasons which we shall explain in a moment, including lack of usable formulas. 
In order to evaluate the inner product (F,cf>j^), one needs the spherical harmonic expansion of F; then 
one uses g(tM)F — g(tJ2i lPi)F to evaluate the inner products. (See the bottom of page 9 of [12].) In 
CMB there is a large region of missing data on the sphere, called the "sky cut" , arisng from interference 
from the brightness of the Milky Way. Thus, finding spherical harmonic coefficients, which depend es- 
sentially on the global behavior of F. is problematic. In contrast, when Mexican needlets are used, one 
can effectively evaluate the inner products if Xjj is somewhat away from the sky cut (on the scale j), 
because of the Gaussian decay of the (fjj at each scale and the lack of oscillation of Mexican needlets. 
(Here we are of course assuming that the approximation (69) is rigorously justified.) 

As to whether effective formulas could someday be found for needlets that could be used directly on the 
sphere, we are pessimistic, for the following (admittedly circumstantial) reasons. If one takes the inverse 
Fourier transform of a cutoff function, which resembles a characteristic function, then one would expect 
an answer which looks like the inverse Fourier transform of a characteristic function, i.e., something that 
looks like the oscillatory function sin(x)/x (of course it must be in the Schwartz space). One would 
expect something similar to happen on the sphere, and in fact, if one takes t = .1 and g(s) to be the 
cutoff function exp(— 1 / [9/ 16 — (s — 5/4) 2 ]) for s e [1/2,2], Maple says that 47r/i t (cos#) is as in Figure 
2, which should be contrasted with the better-behaved function in Figure 1, where we took t = .1 and 
/(s) = se-.) 

Note also that, since a cutoff function is not real analytic, its inverse Fourier transform cannot decay 
exponentially; one would expect something similar on the sphere. Also note that on the real line, if one 
wants a function g such that both g and g are very small outside compact sets, the uncertainty principle 
says that a Gaussian is the best choice; one would expect a similar phenomenon on the sphere. 

Errors of one sort or another being unavoidable, it would be worthwhile to utilize both needlets and 
Mexican needlets in the analysis of CMB, and the results should be compared. We also suggest that a 
"hybrid" approach be attempted, combining the ideas of this section with those of [14], [15]. Instead of 
using a cutoff function, we let f(s) = se~ s , let K t be the kernel of /(t 2 A), and again define (f>jj by (71) 
(for suitable We would then expect (from (69)) that we will be able to evaluate the (F, <j)^i) without 

first finding the spherical harmonic decomposition of F. However, the {4>j,i} are then only a nearly tight 
frame, for the following two reasons: 
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Figure 2: A wavelet on S* 2 obtained from a cutoff function (t 
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a) We have only (3) and (4) instead of point #1 above. 

b) / does not have compact support, so that cubature formulas will not exactly hold. 

However, these issues should lead only to very small errors. For point a), we recall that B a /A a — > 1 
nearly quadratically in (4). We simply need to use dilations by a? instead of 2 J for j sufficiently close to 
1. (As we have said, if a = 2 1//3 , then B a /A a = 1.0000 to four significant digits.) For point b), we note 
that / has exponential decay at infinity, so that it basically has compact support for all practical purposes. 

Let us outline an argument to explain why the {<frj,i}, as in (71), associated to f(s) = se~ s are extremely 
close to being a tight frame (for suitable It would not be at all difficult to make this argument 

rigorous. (In the terminology above, we shall restrict attention to M. — We will pay a price in that 

l(j) will need to increase with j, ever so slightly faster than in [14], [15]. 

Before we begin this argument, it is best to give some elementary estimates that we shall need. First, in 
(1) and (4), with f(t) = te~*, note that c = J °° te~ 2t dt = 1/4. Also te~ 2t is decreasing for t > 1/2, and 
ifM>0, 

On the other hand, if C, b > and g is a continuous function such that \g(x)\ 2 /x is decreasing on [C, oo), 
then by the method of the integral test, 

£ (\g(ba 2 ^\ 2 /a 2 ^-a 2( ^)< / (\g(t)\ 2 /t)dt. 

{j:a 2 U- 1 )>C/b} C 

Applying this to g(s) = se~ s we find that 

{j:a 2 <3-i)>M/6} V 7 

if M> 1/2. 

Seoondly, for j e Z, let j_ = max(— j, 0). Say a > 1, N,l,r > 1; we shall need a very crude estimate 
for the solution m = m(N, r, I) of a 2m l(l + n — 1) = N + rm,-. (Note that the left side increases as m 
increases, with range (0,oo), while the right side is positive and nonincreasing as m increases; so there 
is one and only one solution.) Let V = \ log a (iV/[Z(Z + n — 1)]). If 1(1 + n — 1) < N, then m must be 
positive, and hence equal to p. If 1(1 + n — 1) > N, then m must be nonpositive. In that situation, 
a 2p l(l + n-l) = N < N + rp-, so m > p. On the other hand, let q = \ \og a (rN/l) = \\k>g a (rN) -\og a l], 
so that q_ < \ log a I. Since xy > x + y if x, y > 2, we have 

a 2q [l(l + n- 1)] = rN(l + n- 1) > N + r(l + n-l) > N + r log Q I > N + rq_; 

so to < q. Altogether, we always have 

hog a (N/[l(l + n-l)}) <m(N,r,l) < Uog^rN/l). (73) 

We now present our argument. Since, in (4), c = 1/4, we may choose a > 1 sufficiently close to 1 that 
> ^k^i B a < §t^ 5 where e x « 1. Say again that (I - P )F = F. In place of (70) we have that 
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for some r\\ with \r]i\ < e\, 

1/2 / \ 1/2 



/ \ 1/2 / oo \ / oo 

(^tS i|f|I = ( E i/(« 2JA )i 2 ^> = E n/(« 2jA )^i 
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We repeat that a has now been fixed; we next need to choose N > 1/2 large (how large will be explained 
in the following argument). For j e Z, let j_ = max(— j, 0), r = max(l,2(n + 2) In a), and set fij = 
X[o.Na 2 +ra 2 (j-i)-].f $ij = f — fi,j- Since an I 2 norm is a norm, we can write 

\ 1/2 / \ 1/2 

OO \ / OO \ 

£ ||/(a 2 ^|| 2 = ^ ||/ij(a^'A)F||2 ^ 

/=-oo y y.3=-cx) y 

/ x 1/2 

where E x < (Ejl-co WhA^ ^) F \\l) := .4(F), say. Since f hj (s) vanishes for s > No 2 + ra 2 (j - 1)_, 

/ij (a 2 ^(7+n-l)) = if Z > + r(j - l)./^'" 1 ; let /(j) be the least integer greater than y/N + r{j - l)./^'" 1 . 
Thus 

\ 1/2 / \ 1/2 

OO \ / OO 

£ ||/ij(a«A)F||i = J2 T, X »U)MiA* 2i &)F\ 2 (x2i {j) ,i) 

i=—oo J \ j — — oo i 

Letting K t (x,y) be the kernel of f(t 2 A), we see that this equals 

\ 1/2 / \ 1/2 

OO \ / OO \ 

E E a 2k,),i/(« 2ja )^i 2 (^o),) +£ 2 = E Ei^fe)i 2 

/=— 00 i J \j— — oo i J 

where <j, jti (y) = y/X^K^x^y), and E 2 < [E~-co E, M^lhA^ A)F\ 2 (x 2l{]) ^/ 2 := 
say. If we can show that (for N sufficiently large) that A(F),B(F) < e2||F||, where e 2 << 1, then by 
combining the above facts we find that 



1 + 772 , 
8 logo 

j = -oo 



f\\1= E Ei^fe)i 2 



for some 772 << 1, as we wanted. But 

CO 

A(F) 2 = < £ |/ 2j (a 2j A)| 2 F,F)< e3 ||F|| 2 

J=— oo 

if 

OO 

e 3 = max ^ |/ 2jJ (a 2j 's)| 2 < max ^ |/(a 2j s)| 2 < Ce"^ << 1 

j=— oo a 2 i>Na 2 / s 

if is sufficiently large. (We have used (72) and recalled that N > 1/2.) As for B(F), denote the 
kernel of / 2lj (* 2 A) by J t (x,y) = Jf{y). Then for any x e S n , \ f 2J (a^A)F\ 2 (x)] < \\ J* f 2 \\F\\. For any 
orthogonal transformation T, J t (Tx,Ty) = J t (x,y) for all x, y; accordingly || J^|| 2 is independent of x, 
so it equals 



'n 1 / ||J^||i^(x)=a;- 1 ||/ 2 ,,(a^A)|| 2 , 
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where the last || || 2 denotes Hilbert-Schmidt norm. On the other hand, for any j, J2i ^2i(j).i — J2i ^2i(j),i^ — 
u n , since 1 e Vuj\H. Putting these facts together we see that B(F) 2 < e4\\F\\ 2 , where 



£4 = E \\h,M 2i ^)\\l= E El-M^ + ™-l))| 2 dimW ; 

j — — oo j — — oo l — l 

<cf> £ \hA^Ki + n-i))\ 2 

l—l j — — oo 

oo 

< C E r E |/(a 2 ^(/ + n-l))| 2 . 

1=1 a 2 i>(Na 2 +ra 2 ( 3 -l)-)/[l{l+n-l)] 

In the inner summation, a 2 ^ ^(7 + n — 1) > N + r(j — so j — 1 > m(iV, r, (notation as in (73)). 
Accordingly, in the inner summation, if we set m = m(N, r, I), then 

a 2 ^-^l(l + n - 1) > a 2m l{l + n - 1) = N + rm_. (74) 

Note also that, by (73), N + rm- > N + rc(N, I), where we set 



c( Nl = (1 ° ga 1 l0ga ^ ^ + n ^ ~ N ' 
1 otherwise. 

By (72) and (74), then, e 4 < CJ^Zi l n e' N ' T < N - l \N+r[m{N , r, Since ex V {-\ log a /) < a -(" +2 ) lo s» ' = 

/-™- 2 , we easily find (by (73)) that 

e 4 < Ce' Jv iV™+ 3 « 1, 

as desired. 
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